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ABSTRACT 


Next nearest neighbour interactions are introduced 
for two different lattice models. The canonical partition 
function of the one-dimensional spin 4% ISing model in zero 
freld with nearest (J,) and next nearest (J,) neighbour 
interactions is evaluated. The locus of zeros of this 
Parmtat1 One function is. then. found, for different ratios: of 
J5 to Jy - We also consider a one-dimensional Heisenberg 
model in zero field with nearest and next nearest inter- 
actions, Jy and J. respectively, as before. Repeating the 
work of Niemeijer, a temperature-dependent Hartree-Fock 
approximation is used to obtain a trial Helmholtz free 
energypem site, . whichs-iLs an upper bound. fox the: zrue: tree 
energy per site. The resulting equations can be solved 
at zero temperature for an approximate ground state wave- 
function; the corresponding energy per site is an upper 
bound for the exact ground state energy per site. Several 
cases are looked at, for different relative values of Jy 
and Jo: For J, >0 part of the upper bound for the ground 
state energy per site which was overlooked by Niemeijer 


is obtained. Also, asymptotes to the energy upper bound 


LOX Jy greater than and less than zero are found. 
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Chapters (hs Introdiwcetion 


Much work has been done in the past fifty years on 
the development and solution, both exact and approximate, 
of various lattice model systems in one and higher dimen- 
Sions. Nearest neighbour (nn), next nearest neighbour 
(nnn) and even higher order interactions have been used 
in the study of lattice systems. Models have been in- 
vestigated for both zero and non-zero magnetic field (#). 
Unfortunately, however, most of these model systems have 
not as yet been solved exactly and thus answers to many 
questions concerning them are either inconclusive or un- 
known. 

Ising? Jo yea the mega named after him with nearest 
neighbour interactions and for all H in the one-dimen- 
Sional case in 1925. Onsager's famous solution” for the 
two-dimensional Ising model in zero field with nn inter- 
actions appeared in 1944; however, the three-dimensional 
problem is yet to be solved. The XY model Has been solved 
only in one dimension?’ ie, and that only for nearest 
neighbour interactions with A= 0. No exact solution has 
yet been found for the Heisenberg model in any dimension. 
A happy exception is the spherical model which has been 
solved” in one, two and three dimensions with nn inter- 


actions for all #; however, the spherical model is not 
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Very soauLs VIng.On Warious physical grounds. 

Introducing next nearest neighbour and even higher 
order interactions produces a model which is closer to 
physical reality because there is no physical system 
whose statistical and thermodynamic properties can be 
completely explained by forces of finite magnitude which 
are so short-range that the only non-zero interaction 
energies are for nearest neighbour pairs of spins. Of 
course there are some physical systems which, to a good 
approximation, can be described, over a certain tempera- 
ture range, by a lattice model using only nearest neigh- 
bour interactions. An example of this is the spin % 
Ising model of a binevyealloy. 

The model systems mentioned, TSing, <., ea senberg 
and spherical, all involve a spin-spin bilinear exchange 
interaction. Work has been done on other more general 
Piteract tons such acwthe spin-spin biguadratic exchange 


of is the Dzyaloshinsky interaction® and the 


interact Lon 
dipole-dipole interaction’. It is even possible to post- 
ulate a tensor-tensor type of interaction. All of these 
more general interactions are more physically and experi- 
mentally relevant than spin-spin bilinear exchange inter- 
actions) but, unfortunately, are also theoretically in- 
tractable at the present time. 


In this Thesis we introduce next nearest neighbour 


interactions for two different one-dimensional lattice 


ho 


models. For the spin % Ising model in zero magnetic field 
ther; canonical partition Tunctionsss found™=and its locus 

of zeros is then obtained for various ratios of the nnn to 
nn interaction energies. We also consider the zero field 
spin % Heisenberg model. Repeating the work of Niemei- 
Here we obtain an upper bound for the exact ground state 
energy per site at zero temperature for various different 


values of Jy and J the nearest neighbour and next neigh- 


9! 
bour interaction energies, respectivelys”’ Part of* the 
upper bound obtained by Niemeijer for the antiferromag- 
netic case, namely J,70, is shown to be incorrect. 

It is a well known result that one-dimensional 
SVGeeMeomidving, Short range forces! of finite magnitude 
Cannot exhibit’ phase transitions with’ a. transition” tem- 
bebacune, Tee greater than zero, though one can have 
T 70 with Ong ee nae forces. Introducing nnn interactions 
into systems with long-range forces will alter the transi- 
tion temperature. Bringing in next nearest interactions 
into our two models, even though it produces no change in 
Te Since the forces involved are of short-range, makes 
them more physically relevant and the resulting systems 


possess interesting features which make detailed investi- 


gation worthwhile. 
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Chapter II. One-Dimensional Spin % Ising Model 
1. Preliminaries 


A useful method in discussing phase transitions is 
to investigate the distribution of zeros of the partition 
function. Physical and thermodynamic quantities can be 
determined once we find the locus of zeros and the distri- 


button function on. 6. Suzukaoo 


has derived the asympto- 
tic form of the Helmholtz free energy and magnetization 
for the Ising model. in terms of :the distribution function 
of zeros in the complex fugacity plane. Even earlier, 
Yang and maeaae in discussing the phase transitions of 
lattice systems, proved that phase transitions occur 

when the locus of zeros of the partition function crosses 
the positive real axis of the complex fugacity plane with 
finite density .of .zeros .(density is always) zero on the 
real axis) as N > », where N is the total number of latt- 
ice sites (Yang-Lee Theorem). These are motivations for 
the general .study. of -zeros..oG ppartitron, functions. 

There are two specific reasons for looking at the 
zeros of the particular partition function examined in 
this. whaptes... Facst, Stephenson’? has shown for a cer- 
tain range of nnn to nn interaction energies the existence 


of wm disorder pointy, T for which there is ferromagnetic 
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short-range order below this finite temperature and oscill- 
atory exponential decay in the correlation above it. We 
want to see if the existence of a disorder point is reflec- 
ted in the locus of zeros of the Pave LLOn LunceLon . 

. Since the model we are investigating can be solved 
exactly, one can justifiedly ask why study the zeros at 
all. The answer is that our calculations are being done 
as a trial run. = Zeros, for fixed “ratios “of Jo/Iy have not 
been studied before; previous investigations fixed the 
fac1oO Of Jy and Jo tor ahe Cenmperaturg. iP. Jones*4 pros 
posed a picture of a phase transition using complex tem- 
peratures but the complex fugacity plane is better suited 
SO Lees .Or wweLOs.Ot Daim LOM functions form fixed 
ratios of J5 co Jy. The second reason for studying the 
zeros of eee Simple one-dimensional model is, therefore, 
that it gives fe an appreciation of the mathematical and 
computational intricacies involved in studying zeros from 
the viewpoint of fixed ratios of Jo/5,- Such an, apprec- 
jation is necessary before tackling the same problem for 
more complicated one and two-dimensional systems. 

In this section we consider a spin % Ising linear 
chain with nearest and next nearest neighbour interactions 
present. Montroll?> has obtained the partition function 
of this system for the case where the magnetic field is 


zero by the method of an eigenvalue problem while Obokata 


and Oguchi?® got identical results by generalizing the 
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Bethe approximation. cee rCunc) fhe, partitions tunc— 

tion in the presence of a magnetic field by obtaining the 
eigenvalue equation of the transfer matrix for the system. 
mhe locus of zeros of the partition function for non-zero 
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magnetic field was examined by Katsura and Ohminami' =~. In 


WieterOlbows, alter cbtaining the; partition function by 

the transfer matrix method, we then investigate its locus 
of zeros in the case of zero magnetic field for different 
ratios of the next nearest neighbour to nearest neighbour 


interaction energies. 
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Z2erneros Of PartLltion Function 


Consider an Ising linear chain with N lattice sites, 
having Jy and Jo as the interaction energies between 
nearest neighbour spins and next nearest neighbour spins, 
respectively. Let #H be the magnetic field and yw the 
magnetic moment associated with each spin. It is conven- 


ient to introduce a set of N dichotomic spin variables, 


(oe eee wiere mu Cheats th) lattiseicite, 


cris Que pe: Cyorci@ cube: 


—|, £02 Spin down 


Fo a tpaneroular comtiguratiwon of spins {s}, the 


Hamiltonian is given by, 


N=) N=Z N 
ee rc gees ore) at) ca fey 
i=l i=l i=l 
Introcuce thesnotation, 
Ky = on: ; Ky = ne) 7 obe= Ht (2-2) 
kT kT ae 


where k is Boltzmann's constant and T is the absolute 
Lemperavure, 
The .Canonical partition, function for the system then 


1S > givenyDby., 


OS re \ exp (-H/kT) 
{s} 


aN 


At i N-1 hee 
eu ) exp ( Ai ) iChat lente ee ) ited 2 
o,=-l Oyewl i=l i=l 
N 
pe can (223) 
1=1 


We consider the casé#H = 0. Thus (2-3) becomes, 


i) 1 N=1 
Bike Roy Olea) be) CAPES |) GG. 
o,=-1 Oya i=l 
N-2 
Tae) SAG | (2-4) 
i=l 
| pile 
We now quote Dobson's proof ite te, 
| Z (Ky 1K, ,0) = 224 (Ky 10,K,) (25) 


Introduce a new set of N variables defined by, 


SO): (2-6) 


Ohree telt 8} te We ea NL (2-7) 
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The inverse transformation is, 
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The set of variables fa, },i=1,2,...,N have possible values 
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Transforming the sums over the o's to the corresponding 
sums in terms of the new set of variables ta.t, Our partl= 


tion functrvon (2-4): can be written, 


irae 1 N-1 
Zag (Ky 1 Ky 7 0) ) ), ‘, exp (K, \ OF 
Og-"-lLaj =-l oa. y=-l i=l 
N-2 
HK et One | (2-10) 
i=l] 


The sum over a, can be performed immediately to give a 


0 


result in the dummy variables Ole which can be recognized 
as 22%. 4 (Ko 10K), thus: proving (2-5), The signaticance 
of Dobson si result. Ss. thay ite reduces, the~partitaon func-— 
tion of a linear chain with N spins having nn interactions 


J, and nnn interactions J, and magnetic field H = 0, to 


uy 2 


twice the partition: function offfa, linear chaim/with N-1l 


spins having nn interactions Jo and no nnn interactions 


anidumagnebice field Ht = Ji/u- 


Changing the variables in a result given by nets 


Die oe Une) aes ets (le as) (2-11) 
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where V is the following 2x2 matrix, 
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ak2thi a2 
V= (2-12) 
an ek27Ki 
V is the transfer matrix of the system. Now, uSing a 


theorem from linear algebra, 
Tr(a) =D Fea (2-13) 


MeRaGe, She lhWOMClOen values Ob The strane rer 
Machu, | PRe Choaracteriocic, or secttlar,; equation. of 


Ves, 


eae coeh hye tates eee ene 01 | Lada) 
from which we obtain the two eigenvalues to be, 

ee = e2eosh kK. + (ue See + Spey (2-15) 
ey Wein U2); 

Te RAO Oe aia a) (2-16) 


Equation (2-16) is the partition function we set out to 


evaluate. Its locus of zeros will now be examined. 


Setting (2-16) to zero we obtain, 


{ aathenge try 


o (2-17) 
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1/N-1 i (2k+1)q7 
A ss exp (RS ) 
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The allowed values of k are restricted in (2-18) because 
we Want .only those (N-1)st roots of (-1) which lie on the 
first Riemann sheet, namely between 0 and 2m radians in 
the complex” number plane. Using (2-15) and (2-18), (2-17) 


becomes the following, 


anne Ee + ares = -tan“$,, (2-19) 
saHe 
where, 
o, = pez (kts); k=0,1,...,N-2 (2-20) 


The transcendental equation (2-19) gives the locus of 
Zetes Chane Pod i ton Lunctson 2—16)s, 
For a particular system, Jy and J. are fixed real 


numbers. Therefore we can define, 


A =J > A real; —-e<)\<to (2-21) 


Note that if we replace TT by “Jas (2-19) remains un- 
changed." "This’means~that. the roots"of*’ (2-19) are not 
affected by a change in sign of the nearest neighbour 


interaction energy, Jy- Therefore, without loss of 
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generality, we can set Jy = FL. com: (2-2) 5 

oe SASH: (2-22) 
and, 

A = K5/Ky 7 | (2-23) 
Oi, 

K, = AK) = /kT (2-24) 


Using (2-24), (2-19) becomes, 


-tan dy (2-25) 


Werdetine the “fugacity, 7. 2, by, 


—~2Ky 


z=e (2=26) 


Strictlywspecking, Zz 1S? not,a wrugacity but 1t' as the quan- 
tity analogous to the fugacity in the grand canonical 


ensemble. In terms of z, equation (2-25) becomes, 


cos“$, = (az) ; zt0, 27°41 (2207) 
a7 i Zee) 
The roots of (2-27) are the zeros of the partition func 
t20n (2-16). the amaginary part of, (2-27) igives the locus 
of zeros and the real part gives a subsidiary condition 
that the zeros must satisfy. Numerical results of our 


calculations are given in the next section. 
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S.e Numer veal Results 


To facilitate computations we introduce polar coor- 


dinates (o0,¢) such that, 


Z= pe (2-28) 


The advantage of using polar coordinates in discussing 
(2227) is tthat, -y So, doing, we canbe sure thatewe: will 
obtain the roots on the first Riemann sheet and thus avoid 
the complications that come with multi-valuedness. Taking 


the: complex conjugate of (2-27) we obtain, 


Hea ze Oy x7 Aty (2=29) 


cos“$,, = 
42% (1—-2*74) 


where..z*,as* the, complex, conjugate’ of: 7. . Subtracting 
(2-27) andy 2529) 7 “then, bringing) the. difference) to, a 
common denominator and using the four assumptions z+0, 


Sat eto, 2% 744) wer yeinds that: 


ee) Sane - {sin (2041) 6 + 2o0sin2do 


+ p-sin(2-1)o}= 0; 2#0,27481 (2-30) 


Equation (2-30) gives the locus of zeros in the complex 
fugacity plane. 


Now, in equation (2-27) we multiply through by 


2) 


4(1-z); we can do this since we have assumed ye Mit 
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Equating real and imaginary parts we obtain the following 


two equations, 


4cos*$,, = (ge eoyooss ico. 1-p~*cos21$40 


to" cos2\A¢ 


ety pez 
and, 


2r 


4cos*$,, = fie Saeine a0) sin2Ao40,27*41 (2-32) 


pe 


BOouations (2-31) and. (2-32) are two additional conditions 
that tthe mootsito® (2-2) )*mist: satisfy7r Of “the three *equa-— 
tions (2-30). 4(2-—31)mgand: (2-—32)2 only twoi-are andependent, 
as (2-30) can be obtained from the other two by the elim- 
ination of "4cos*$,,". By Padding) (2=31)) sand. v@ =32) an 
equation can be found which is directly related to the 
subsidiary condition mentioned at the end of Section 2. 

We wail refer to )(2-31)-and (2-32) as two subsidiary con- 
ditions that the roots -of (2-27) must satisfy, even though 
this terminology is not strictly valid. 

We are interested in the case N + », Then, dy is 
Continuous oan the anterval O<o, ST and cos“$, is continuous 
in the interval {0,1}, where the curly brackets indicate 
a closed interval, namely one in which the endpoints are 
gncluded. ‘Thus, in order for: the subsidiary. conditions 


(2=3)) and (2-32)\ to be satisfied, the right hand sides 


' 3 9A3 rita ido Ye 


of the two equations must be in the interval {0,4}. 

We have solved (2-30) numerically for various values 
of A, subject to the constraints imposed by (2-31) and 
(2-32). The results for the locus of zeros of the parti- 
tion function (2-16) in the complex z plane are given in 
Figures 1,2,3,4,>5 and 6.. 

ti oy (2-27) we make fhe transtormativoen Aes =), 
z+ 1/z', we obtain the same equation in the primed var- 
iables as we had in the unprimed variables. This means 


Eee Neicer OL roots Of (2-2 /)e corny = AG eS. vans 


Ghen, the corresponding Set or Toots Of (2-2/7) Tor*x% =-2>5 
fo 17 25) venue, co ,find theryroots for’) ="-\5) alt we 
need to do is a complex inversion of the roots for hd = do. 


Because of this symmetry in (2-27), we have not separately 
considered any negative values of X; positive values only 


have been investigated. 


” 


a, = 7 
; Oa avy 03 


esufsy enoigne wt y Ph 2 te 
bas (£2<0) qi Beangal 


cn 
~igweq ai? 34 *6ene GA 2768s 


] 


nh qeviy oxs aonly s <c ogoee 


+ 


aa 
oO 7 


AO Taeyes inst oe Lites sw (vt > te 

7 . - 

—tea ft 34% thLonei~ipe:amae art ulnddo au 4 z 
mi y 24ts Fp aang age al Bat ae 


sa : '"3<°) Yo! adept go doe ant) t 


. 
, 


(v4 ; es & ie, 28 erkbooges1700 ond, f 


42.0% of , endT 


ef Fi. 
on" i- 


rit es { of (1M, a) el ob 


| 2h YSOGmeye elias 


[ow evéernpad Yre 


_ 
,e@7 Br > severtek steeped: . : 
7 


Pvoures) W275 ,4,0e¢0nd 6: LOCUS: Of “zeros of the 
Dartitaonm Luncti1on (2-26) in the complex 
fugacity plane for different values of 
we bor all graphs horizontal exis gives 
realspart of z and vertical, amismsthe 
imaginary part. Scale for two axes is 
the same for all twenty-two graphs. If 
ea particular~branch of Locus” of zeros 
goes beyond the scale we have used, the 
point it approaches is given in terms of 


polar coordinates (p,6) . 
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0,p 2216 : 
A= aT pike. 
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A. Dascussit6on Of Numerical Results 


A. Properties Of Solutions - The General Case 


We now work with equation (2-27) to obtain some exact 


results by analytical means. First we investigate solu- 
LLOns, fOr oO = 102 (2-27) becomes, 
2 
cos hy = SG mergoto, £1 | (meen 
2A 
AO (e'pe #) 


We want those real positive values of 9 for which 
0<cos“$, <1. USing (2-33) we find ‘that therevare no 
values of 0 which satisfy both of the inequalities simul- 
taneously. This means that there are no solutions of 


(2-27) [or oO. = 0%, 


FOr 0 = 1; 
Z Z 
cos oy = (sO) es > pt0 
~4p (1-p2* e714") 
: pa eet ay. (2-34) 


Using the continuity, of cos*$, in the closed interval 
{0,1} in conjunction with (2-34) we can examine the three 
cases 2). an even or odd integer and 2r\¢integer. We. get 
the following as solutions for 9 = T: 


a). l<p<~™, for 2X = even integer (A=1,2,3,...). 
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b). p = 1, otherwise. 
Our numerical results, as given an Figures 1 through 6, 
bear this out. 

Next we look at the situation when ¢ = 27. Equation 


(2-27) Gives’, 


cos“, = Ele Meo Oupg ea er (2245) 


Ney Gl Syeieg es A 


Again we want those real positive values of o for which 
0<cos“$, <1. For 4X) an even integer and 4\#integer, (2-35) 
has no valid solutions; however, for 4X an odd integer 
bel /4+3 (4, S Aspaeed) pythere eas o interval for which the 
requirement on cos“ by is met. The numerical values of 
the first four intervals are the following: 

ewer 01.92.6333 0 50419" 47.05 i cebor é) <2 .1/ 4% 


Db) aor 0238 L849 ~ifow Agr 8/4. 


@)istoe0ie4 G8S00:., «for x 5 Ai 

Gh. G2 Or. 5126 Fl, GEGIn wcsie] AA: 
Intervals BLO yi = 9764) We ay 1G) DeSrWwOLKed: Cutealso. 
Thus, solutions for 9 = 27 exist for odd positive quarter 
integer values of \ and the solutions are either finite 
Quisel (A) vomeseni-imtinite ((\=3/74,5/4,...) intervals 
along the positive real fugacity axis. These solutions 
are shown for the appropriate cases in Figures 1-6. 


We can-also examine solutions of (2-27) for op = 1, 


for which we have, 
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cos“, = ie ; eli aa (2-36) 


Sampisrying (2-36), then separating real and tmaginary 


parts we obtain, 


l+coso 


= cos”, ; Sep ee ey (2-37) 
4 » 


and, 


(l+cos¢) sin2)d6 
4(1-cos2)\$) 


0 en, tome cae (238) 


We want only those roots of (2-38) which satisfy the 
constraint imposed by (2-37). For the first Riemann sheet 
we consider only those values of ¢ in the closed interval 
PO i) ae BOG levalues Of On 1 wtntowinternval, cos}, 
given by #(.2-37) Ssatisfiies 0<cos“$, <1. Thus the noots of 
(2559) are the colutions OL (252 We tOr nO ule We find 
Eoe. Solutions, Lorso = 1 sto be. the Col towing: 

Spd eet Or oan teqer a | ly easy wen)). 

b) . o-= oy : m=O, eee es Pa 7 aeOr ay +, integer 

enon Olen Od. 

Che. Hao petom AN = Oddhinmteger = 1/4, 3/4 5/744 .3\5.) 

where we have defined, 
greatest integer contained in 


Baise See Or se 0 


Lp 
0, otherwise (2-39) 


Tie first. and. third cCypes: of “solutions for 0 = 1 ere 
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consistent with the solutions we have already found for 
¢ = 7 and ¢ = 27, respectively. As an example let us 
Sonsacdere\t-391525.. ,Tlwossolutionsson (2-27)) ané ove) = 


Clit bland &toy0) f= (1527). deSimee [4)] = 5, the allowed 


ae 
2) 


an integer, so m = 5 must also be excluded. Thus the 


Gade valuessoremmanesina= bF3 [5 segetoweves).{ for mn= 15% = 2, 
Geher; two solutions ofe(2-27)mitorsorS heer i= 19825 are, 
Hromim =f) fand73,*40,¢). t= (1,60), (1,2m). These four 
Ssobutionspomes icanébetseenuiny the tgrapWifor Ni=ube25. in 
Figure 5. 

Nest; iwer “look Sir Slange *oNsclutvons of (2-2). "Kor 


solutions where z is large in magnitude, . 


cos” $y, aL aerate es (2-40) 
4z 


Since cos” $y must satisfy 0<cos“$, <1, this means that 
Mange. e soluterons are possibleioniviast Qaairds thatrais 
Ns 3 Le gwa tingi iieall Vandi Imagineany parts, ini ,(2-40.) ; 


oy 


and, 


Saree) yo aes pao, ee (2-42) 
at) ee 


Solving (2-42) subject to the constraint imposed by (2-41) 


anc Weng the fact that 6, is restricted to the first (Rie- 


DET | 


mann sheet, 0<¢<21, we arrive at the result that (2-27) 
has* large*ovsolutions; namely o sPer,oforad = n71/(2iA=1) 9 
where n = Teoe5pets a=) ast 4,-2] PSR Odduand. Nis. ay.3', 
5,...,[4r-2] -1 sh f4r-2] is even. From this it is easy 
to show that the number of large op solutions in the 
closed interval {0,27} is [2a-%] if [4-2] is oda ana 
[2-1] oe [4-2 is even. Thus there are large p» sol- 
HtP#onSssoOLthe-27jutoreall values ofypu satisfying 2a—s21, 
namely for all 4123/4; our numerical results, Figures 
136 thea. thi shout. ee ete usilooksatacthe, case ) = 6. aS. an 
example. [4a-2] isMevernusoasunere. shoudd be [2a] = 5 
large p solutions, and there are. They should occur at 
Oeew ny S44 ne=i1el 8h oy 70 Wohath us St 6Perw/55 3775 ha 175, 
Ong ee Bron thergraph for) =spnmehigune «6 47 darge\, o=sol- 
utions do exist for these five values of 96. 

We note another interesting feature of (2-27). If 
(ond) =88tp4 ,7-oqavisma solution, thentt(o.,6)o=9c(o43, 746, ) 
iS"olscoza Soluiscrewhen O28 = n72i-sree) 0341, at. o THIS 
symmetry is reflected in the appropriate plots in Figures 
1-6. For no other values of i} other than integer and 
half-integer values are the loci of zeros symmetric about 
gd= 1. 

Finally, we give a qualitative description of the 
behaviour of the locus of zeros for increasing ); this is 
in reference ‘to ‘the graphs in) Figures 1-6. For » small 
and increasing the locus expands out from the point z = 


~l1. The branch in the upper half of the complex z plane 
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bends toward the > = 0 axis while the branch in the lower 
half bends counterclockwise and approaches the positive 
real axis (¢@ = 217); it reaches it when Meee. 2 dae USO 

TOL thus value Of A’ , there 1Ssua finite. Solution interval 
On the positive real fugacity axis. For A>0.25, the Locus 
assumes a roughly circular form pu a bulge which is the 
departure from circularity being below the real axis. At 
hee Cle Blocis 1OOKS) like at circular rarce, —oranc0.5, 
the bulge appears again, this time above the real axis 

and travels clockwise with increasing \}. At dX = 0.75, our 
first large p solution appears at ¢ = 2n from the second 
Riemann sheet. It rotates clockwise and reaches the neg- 
ative real axis (¢ = 7) when } = 1. For A>0.75, the bulge 
changes direction and again approaches ¢6= 7; at » =1, 
Chis eoranch of the. locus again looks: like avjcarecular arc. 
For \}>1, the bulge moves below the real axis and the 

first large 9 solution continues: to rotate clockwise. At 
4 = 1.25, the second large po solution appears at $9 = 27 
and begins to travel clockwise as } increases. At A = 
1.5, the bulge again crosses the real z axis: The third 
large op solution enters the first Riemann sheet at $ = 27 
when \ = 1.75 and also moves clockwise for increasing i. 
The bulge crosses the real axis at i} =2 again, at which 
point the second large p solution has reached ¢ = 7. As 

A increases, we see that we get additional large o sol- 
utions appearing at » = 27 each time we reach an odd 


quarter-integer value of ); these large op solutions then 
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rotate clockwise with further increase in A. One large p 
solution crosses the negative real axis whenever A takes 
on an integer value. More and more branches of the locus 
of zeros appear in the first Riemann sheet as A gets 


larger. 
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B. Exactly Soluble Cases 


For certain values of A equation (2-27) can be solved 


exactly by analytic means. For ‘S24; 
2 2 . 
Costs nee 3 ZO ye Ze L (2-43) 
4G ( b=) 


This can be expressed as a quadratic equation whose sol- 


HELO 1S), 


Z = L = + — + iB (2-44a) 
2 2(1+4cos o)) 


where, 


B= 2cos$, (2-cos*$,) 7 (2~-44b) 
1+4cos“$,, 
We introduce coordinates (r,s) defined by, 
By =a thas 7 (2-45) 
(2-44a) can then be put in the form, 
CL Ie arte tee | (2-46) 


(2-46). :is.<the equation of a circle with- center at (r,s) = 
(=17 3,0) and NMaving radius, 2/3. “Ehus, the locus of zeros 
in the complex fugacity plane for A = 4% 1s a circulandre 


given by that part of the circle defined by (2-46) for 


which 0<cos“$, <1. The endpoints of the arc are at (r,s) = 


30 


(1/5,42/5) and the angle . = einr (27 /5) ae 162 .O°Nisee 


PLogurets)< 


Next! wer consider )> = tor which, 
Z 2 
Csi, giewiltz Fi maO yp zamees) (2-47) 
AZ, (L=Z) 


Expressing this as a quadratic equation we obtain as the 


solutazony, 


2 
(Cy + Co )/C. * O0<4cos oO, <3-2V2 
—= ~ as *5 . pe ea 2 —_ 
gecinie, + 1(-C,) \/cl = 3-2v2<Acos 9 <4) (2-484) 


-l1-/2 ; 4cos”$, = 3-272 


where, 


A ooh 5 
Cy = Acos o. it (2-48b) 
a i 2 Doon 2 a 
Cy : (1-4cos 1) L6COS dy (2526¢) 
C, = Bcos- 6 (2-48d) 

3 k 


For 0<4cos$, <3-2¥2, the locus of zeros given by (2-48a) 
is a straight line on the negative real z axis such that 
OO ee Os 3-2/2<4c08*$, <4, we use (2-45) to express 


(2-48a) as, 


(tee Di es, 2 (2-49) 


(2-49) is the equation of a circle with center at (r,s) = 
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(i Qetaniera unc. Bedius w./2. Thug, tre locus of zeros in 
the Complex fugacity plane for *\-=' 4 us''a straight, line 
On -the: negative: real axis suchPinae—-<1<-l and a..circ- 
lias arc given by that part of thevcirele defined, by 
(2-49) for ‘which 3-2/2<4cos$, <4. The endpoints of the 
are vare located at (r,s) = (372 oy 778) and the angle 

bo = Se ee = 41 .4° (see Figure 4). 


It is also possible to obtain the exact solution of 


(2-27) for X = 1/4,3/2 and 2. However, these three cases 


were not done as they involved the algebraic solution of 


quartic equations. 
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C. Other Special Cases 


As mentioned in Section 1 on Preliminaries, 
Stephenson*? has shown for ~-4<)<0 the existence of a 
disorder point, Ty: If this existence did peculiar 
things tothe, locity of Zeros ofthe par tltLlon function 
(2-16) for -%<)d\<0, then these peculiarities would be 
rerlected in the loci of zeros: for 0<)\<s, since one set 
of loci can be obtained from the other by a simple com- 
plex inversion in the complex fugacity plane. This symm- 
etry of equation (2-27) has been discussed in Section 3 
on Numerical Results. However, the loci for 0<)\<% are 
not unusual in comparison with those for \2>%, as can be 
seen in Figures 1-6. Indeed, as we have described in 
Subsection A of this section, there is a smooth and pro- 
gressive evolution of the locus of zeros with increasing 
AX. We thus see that the existence of a disorder point 


does not appear to have any unusual effects on the locus 


Of Zeros. 

The case \ = -% has a special Significance. In the 
previous subsection, we solved i = % aera ane showed 
the Locus, of zeros to be a circular arc; thus; X% ='-% 
also gives a circular arc. We have J, = —5|J,|- Using 


2 


Dobson's result, equation (2-5), the variables for the 
corresponding nn chain with interaction J(= Jo) and mag- 
netic field H(= J/u) now satisfy aT SNA ae Tivos. Che 


Gonditson tor the critical magnetic field of an antiferro- 
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magnetic chain’, 


Finally, A = 0.25 also seems to be a special case 
due to the presence of the finite solution interval at 
¢ = 2m7(see Figure 2). Since large op solutions are not 
possible until A>3/4, the reason for the appearance of 
this interval is that it arises Gectkhe of the mathe- 
Matzes°of the problem?-“The’ interval has’ no°physical 
SvGna fr teance since o-="27 rs" nSet physical." "For thrs 
same reason solutions at ¢ = 2n for all are non- 


physical; they come purely from the mathematics. 
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Chapter III. One-Dimensional Spin % Heisenberg Model 
1. Preliminaries 


The Heisenberg model, of which the Ising and XY 
models are special cases, is an interesting theoretical 
model to study because it involves all three Spin com- 
ponents. For this reason it is also mathematically more 
complex. For an infinite linear chain of spins inter= 
acting with nearest neighbours only via a Heisenberg 
interaction, the exact energy eigenstates were found by 


Beene. and eaalane ne: later obtained the exact ground 


state energy for the antiferromagnetic case. Orbach** 
derived the ground state wavefunctions for the linear 
antiferromagnetic Heisenberg chain. The exact partition 
function and wavefunctions for non-zero temperature, 
however, are not yet known. 

If next nearest neighbour interactions are also 
taken into account, the rigorous ground state energy has 
not yet been obtained by analytical means. Various 
approximate methods have been used. In the classical 
approximation spin vectors in the Hamiltonian are re- 
placed by classical vectors; Kaplan“~ is one of those 
who have used this classical picture to examine the 


ground state and its properties. In this approximation 


one expresses the energy for an N spin linear chain as, 


N 2 N 2 
Me) 2 Ose. + J, ) S*cos6. 
be ah ae al 
1=1 i=l 


EN je? (371-1) 
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where S is the length of the classical spin vector, 6 


| 1,i4+1 
the angle between spins at sites i and itl, O. $42 the 
fie 
angle between spins at sites i and i+2 and Ji 155 are the 


nearest and next nearest neighbour interaction energies, 


respectively. The ground state energy per spin then is, 
Zi z 
Si) CA oave ig emacgaee 0 and |g, 1<4lg,| 
85, 
SS Oieacd ai Ca) 
S* (43, + Jo) , otherwise 


where the "+" sign regers to ferromagnetism (J, <0) and 


179) - For J5>0 and 


1184155 | the classical spin chain has a spiral structure 


the "—" sign to antiferromagnetism (J 
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The ground state energy has a kink at Jo = 4S 4 


5 for a “tLixed value of Ji- 


Yoshida and Miwa? and recently Ono” have employed 


when plotted as a function of J 


spin-wave approximations to discuss the problem of ob- 
taining the ground state energy with both nn and nnn int- 
eractions present. Majumdar and Ghosh’ determined some 
properties of finite linear Heisenberg chains and from 
these derived upper and lower bounds for the ground state 
energy per spin as the number of spins in the chain tended 
to infinity. Variational methods have given good agree- 
ment with the known exact results of Bethe for the case 


. : ; ee) 
of nn interactions only, so Niemeijer has made use of 


a variational technique for the case where nnn interac- 
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tions are also present. 

In this chapter we check Niemeijer's calculations 
on the one-dimensional spin % Heisenberg model in zero 
field having both nn and nnn interactions present. An 
“Upper sound “fon, the exactsqrounds state energy. per site 
at zero temperature is obtained for various different 
values of the two interaction energies. Part of the 
upper bound derived by Niemeijer for the antiferromagnetic 
nearest neighbour case is shown to be wrong and the cor- 
rect upper bound, which is an improvement on his results, 
Pesrounc.) Letting a) be the ratio7of nnn ton imterac— 
tion strengths, we also obtain the asymptote for the 
energy upper bound for ferromagnetic nearest neighbour 
interactions having large negative a and for antiferro- 
magnetic nearest neighbour interactions having large 


positive a. 
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We, begin, oun, review of Miemeijer's work? by consid- 
ering a one-dimensional chain with N lattice sites having 
Jy and Jo as the nearest neighbour .and next nearest neigh- 
bour interaction energies, respectively. In zero magnetic 
field and for Heisenberg type interactions the Hamiltonian 
for the system is, 

N 
S5°S441 + Jp 1 85°S54o eee 
where on is the i component of the spin operator at the 
j'th site and where cyclic boundary conditions have been 


and §$ = 9S... -For™che ‘quantum 


taken, namely S ~N+2 roy 


oy et 


mechanical spin 4% case, the S's are the Pauli spin mat- 


rices. Even with Jo = Othe souls fiLcule problem. for 
only the ground state and a few excited states are 
known g 14. For, our .problem we use ja, temperature, depen- 


dent Hartree-Fock approximation and we express the Hamil- 


tonian (3-1) in terms of Fermi operators, which can be 


done because we are considering the spin ¥% case. 


First we perform the transformation, 


oa (a + a.) (3-24) 
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in terms of which (3-1) becomes, 
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where the symbol "t" after an operator refers to the 
Hermitian conjugate of that operator. For the ferro- 
magnetic ground state we can take either all spins up 

Ox all spins down: ati does: not really matter which, as 
long as consistency is retained throughout. In this 
chapter we will consider the ground state to.be all spins 
down. The physical significance of a, is sthet 16 isa 
Creation Operator. JApplied to the j°th site in its 
ground state it creates an up spin, which is an excitation 
ofethe ground stete, “Similarly. a is an annihilation 
operator because, applied to an up spin at the j'th site 
(excited state), it annihilates the up spin and leaves 


a down spin (ground state).. The expectation value “of the 
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number operater asa, gives the occupation number of the 
excited state at the j'th site. 


We make the further transformation, 


i ee! $h 
a, = exp{-ri ) cc, te 
k=1 
3-4 
f a (3-4) 
a, = oc expiqi )se.ceh 
k=] 
Equation’ '(3-3)°1is then written’ as 
N +: +; 
= 1. 
) ((yo3e, 44 + ee pis Can) 
{=e 
a - i e 
+ Jaen ee Cadac ate: ieee 23200420 
= (ot ele, = J5 (cle 2 te aC s)) 
Ley 5 2a kee qa tee mks 
+ 4+ 4S oj (35-5) 
where "h.c." means the Hermitian conjugate of the pre- 
ceding expression. The c's and eas are Fermi operators 
Since they satisfy, 
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anticommutation relations at like sites and Bose comm- 
Utation relations at different sitese i By means of the 
transformation (3-4) we have expressed our Hamiltonian, 
equation (3-5), in terms of a set of pure Fermi operators, 
the c's, whereas before in equation (3-3) it was given 
inecerms ,of the. aS >a -SeCtiiOE Paulions. Also, since 


asa, = ile , from (3-4), therefore the expectation value 


j 
of the number operator ae. is the occupation number of 
the excited state at the j'th site. 
Define the real scalar quantity o, by; 


) = Teo | ama at Ate Pray S| (37) 


Next we make a Fourier decomposition, 
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eee a: explio, jin, 
k=1 
(3-52) 

ee +e | 
eletanN gi) exp{-i¢, jin, 
k=1 


The n‘'s and Ave are a new set of Fermi: operators since, 


vie 
LM My uh = Sh, ac 
(3-9) 
{np poe ae Osi ne } 
KI": kK! 


It is easy to show that, 


f 
+. 


- 


~ 


| Boa as 


- 
_ = 


oo 


teoreds , 1b-E) mot +; , 
: - 

ae en 
} #t .5.5 iotaxegqe tedmur 


> 


> ae 


- 
ad 


rt sty te etst6. Py: 


et be 
— a 


ti up zelsoe Jems oat 6 


f palend 
bone iP 
~~ 
aad 
¢; ue _ 
Ore 2 a } 
7 : 
yt 
{ 
7 
— ry 
fans worls 


PE LE Aya te oy NN (S10) 


By generalizing the orthogonality relation (3-10), 
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exp{io,3} = NA(¢,) (3-11) 


where, 
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Ro) = | | (3202) 


0, otherwise 


A($,) is a real scalar function which is even in by: Our 


Hamiltonian (3-5) can now be written as, 
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(3-5) "1sethat, py se doing; we! have*créated=a-set-oE7en 


No 


(GimE4 We nay Sy both, oe 


’ 


» (i) nedintes vaate 


ie 


~~ 


_ 


) 

0 pie eo 

= ete 
= scr 


t j 7 = AL + ra ‘ f 
| Fr it) 
| on 
fe 3.0 
> 
F Pa thi Bae lus e wed ; 
Rel) 
v - wont neti (c=) niatioos Pits 
Sin 
a) 
4 
“ - 
sf eh | i = " 
ee 
i : : e 
m= a a c + F; 
ab ' A 
a tf 
: ‘{ 
Sat’ { |Saoaeh 4 
PS 7 r 
( =i ; 4“? = s + a? a 3 7 


pe 3 


O12 WR ow 


shevrths We oninges Tseoke 
J - : 


spate en ed , Sartd & 


43 


spin waves whose linear momentum is given by dy in equa- 
taon. (32/7), 20k" titselicbeing®theawave: wector ousince, tit 
N> oo , the momentum varies from zero to some finite max- 
imum value, therefore the corresponding wavelengths of 
the spin waves will go from some finite minimum value to 
infinity. Ordinary spin waves obey Bose statistics but 
the spin waves feed here differ in that they follow Fermi 


statistics, equation (3-9). 


We introduce a trial Hamiltonian, Ho, given by, 


" + 
og) dee) Oy (3-14) 
k=1 
where e(k) is the one-particle energy. The trial Helm- 
holtz free energy, Foy is defined by, 
Fy = <€H>o - TSo (3-15) 


where (H),) is the canonical expectation value of H, 


-BH 
Tr(e ot (3-16) 
Tr (e ©) 


{H)o = 
with, 
eee te 
8 = qT (3-17) 
where k is Boltzmann's constant and T is the absolute 


temperature. SS, is the entropy of the system described 


by the trial Hamiltonian (3-14). If the exact Helmholtz 
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free energy is F, then Bogoliubov's inequality7® states 


that the trial free energy, F is an upper bound to the 


5 aoe 


exact free energy, F, i.e. 
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We can minimize this upper bound by varying the one- 
particle energies, e(k). 
The expectation value of the number operator nny 


is given by, 
n = antn Yo = : (3-198) 
k kk 1408 () 
From: (3-3) (it isweasy, to.show bywusing (3210) .that, 


kk 


N 
y auc =Jni (3-19b) 
j=1 k=1 


The Ws, unlike: the a‘s. and the c's, are not associated 
with individual sites. We interpret n, as the occupation 


number of the spin wave having wave vector "k". By (3-19a), 
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We now employ (3-21) to evaluate <H)») defined by equation 


(3-16). We obtain, 
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where, 
€9 (k) = J, (coso, -1) + J, (cos2¢, -1) (a23)) 


With N an even integer we assume that, 


eS e Rye e(S Pa iy (3-24a) 
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Thus, we have assumed that the one-particle energies, 
e(k), are symmetric with respect to the point oy =v. A 
little later we will show that this assumption is con- 
sistent with currresults. 


The assumptions (3-24) imply that, 


N 
) neeind 0 = ) n, Sin2o, (3-25) 


Equation (3-22) then reduces to, 
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where Ey is the trial energy per site. Now, the entropy 
per.site of the, system of nonanteracting fermions des- 


Criped by (2-14)'5So,,h1ls Given by, 
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where Q,, is the canonical partition function for the sys- 


N 
tem described by (3-14), 


Oy = een 05 | (3-28) 


and E is the expectation value for the energy of the same 


system, 


3 
= Fe alery (Ho) 9 (3-29) 


Combining these results we find, 
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We rewrite (3-15) as, 
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ee jee USK: (3=31)) 


where Fy is the trial free energy per site. By Bogoliu- 


bov' s inequality, Fi LSan, Upper bound to dF) the exact 


free energy per Site, i.e. 
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We see that S is 7 times the magnetic moment per spin in 
thewz direchion.. in the. limit ,aswN +.©, the summations 
in equations (3-34) become integrations. From (3-33) it 
is readily seen that our result fe consistent with the 
assumption SS ee im Ore? = E(d =o, 1) 6 Rovuation (3-33) Le 
the fundamental non-linear integral equation for e(k)-in 


this approximation. 
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5. Solucron For iT = 0 
For T = 0 we have the simplification, 


Oncat e¢k)>0 
ny. = (3-35) 
IF 26 ets) <0 


This allows us: to obtain the exact) results: (3-32) re- 


duces to, 


na > qt — 

Bl Set, Cele), 
The trial energy per site, Eye thus becomes an upper 
bound for the exact ground state energy per site, Eo. 


We are interested in the case N + ~, for which the 
summations in formulas (3-34) become integrations. As 
Ni Seco, dy > @¢ where ¢? iS a continuous variable in the 
closed interval {0,27}. When we integrate, we integrate 
Over ma. Certain interval Of ¢, Canlliic 1s We choose 1 to 
be {x,y} and I lies in {0,n}. Since we assumed e(k) 
symmetric with respect to the point >, = T, equations 
(3-24), therefore I must also be symmetric with respect 


to oy = (see Pagure /)2 Bquationc (3-34) "thus become, 


S = [as - = (y-x) - ¥ (3-37a) 
a 


me) 
II 


2 | cosdas = 2(siny = sinx) (3-37b) 
we 


o is (ce 
4 qu iiés ; 
ey , 
-. i 
( “fy 
i 
ft 
. 
\ iz f SW 
1 = } 
t ee 
rf 5 
’ a ’ ‘ 
' 
i-i] ixut= 


stipes: Jone arta wtnddo cit aii, a 
| I on 


a Zs te 
sy 


{J ste ee (prea oe, 


i eiete boggy daeme, emg = 
) ot) at eenyedal ose ow " 
| a 
i (he iformo? af aaoldenam 
SY. 
rae ‘. <ppes hwy > + a? = 7 
a> 
r Ww 7 | + Nq OG} lavrejal = Lo 
o leven? nie eae 
i 
fae ro ry pL! f Sem iyvixl: ex 7 
. an 
97 IH4Geag ria ry oi wr ia a . 
‘ sauna 7 wr rnegets G8 £) 


an ak eT. | at) YT = ee 
i 7 


! 
f 
°C 
ee) 
i 
ee 
t—— 
4 
o 


= wie) s = phgaun |S =s. 
_ t 


50 
0 = 2| cos2Hag = Sin2y)= ‘SanJx aoc) 
ai 


We must determine the interval I such that when equations 
(oj. are suostitutedt into (G=35) fone nas that. ¢.(k) <0 
when 1. f-and)e(k) +0. when o,€ I. 

There are four cases to considers: (1) J,<0, JO; 
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Our definition of a differs slightly “from-that used) by 


Niemeijer. We shall take Jy = +] when J,>0, and Jy al 


when J,<0. Thee 1s clearly mot.a lamatation. since a) is 


unrestricted in value, namely ~»<a<+o. We can then write 


(3=33:) "as, 
2a 
we ee, 1 2 We 2107) is 
e(k) = (2S+Q)cos Oy + Jy L Fr iger ai COSo) 
25, (1+20)S 
Senior vee enn : (3-39) 
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Equation (3-39) is what we now solve, for different values 


OL 1 A= +1), and for the whole range of a(--<a<+0). 
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ALN FerromagneticuwNN Interactions 


Nearest neighbour interactions are ferromagnetic for 
J,<0 so we take Jy = -] and therefore a = “Jo. We want 


the solution of (3-39) for Jy Snap coc choo (i, Es Oe 
namely X= y, 1s the solution then e(k) must be positive 
everywnere and. from (3-37) 2s = steph =)O = 0.  ihauatien 


(3-39) becomes, 
e(k) = -2acos*$, - cosh, + (1420) (3-40a) 
MoEGD tACEOrIzes as, 


e(k) = ~2a (cosd, = 1) (cos$, te ek. 5a) (3-40b) 


Equations (3-40) are plotted for different ranges of a in 


Pioure 8. For o2-i,ec(k)>0. Por o=-), however, otk) <0 


ae 


over part of the range of d) which means I = 0 is no 
longer the appropriate solution interval. Thus, I = 0, 
hanely = = y, as’ the solution. of (3-29); for a>-3 when 
Jy = =. 


The first factor, an: (03-400) Siacmtwo an zeros at 
$y, =e eis. URLS! NORstNe (oO Tome Or Us= 59) Toren. —o 
because then the second factor in (3-40b) contributes two 
Aa Re Tonal ireat zeros to. ¢(k) “inetne intrerval, 10/27}. 


Cartes, 29 lane en <1 (i <0 nO) el = Ss. the solurion’. 
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Thus, S = <1/2, & =KO_= 0. AS indicated before, S is 7 
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times the magnetic moment per spin in. the z direction. 
In the completely ordered state, the magnetic moment per 
spin in the z direction is -k, since we have taken the 
ferromagnetic ground state to be all spins down. Since 
S = -7/2, this means that this state is totally ordered. 
This state, then, is actually the true ground state. 


(3-36) thus becomes an equality, 


(o- 4 EF) 


The single-particle energy is given by equation (3-40a). 


The exact ground state energy per site thus is, 


Ee EQ = = (1+a) (3-42) 
(25 Jy = =], ped Oda 0, J5>0): Por —750<0) the solution 
is still given by I = 0; the single-particle energy and 


trial energy per site are given by (3-40a) and (3-42), 
respectively. Ey in equation (3-42), however, is now an 
upper bound for the exact ground state energy per site 
rather than the exact energy itself. . 

For a<-k part of the elementary excitation spectrum 
becomes negative and I is now the interval {x,y}, as shown 
in Figure 9. We have restricted ourselves to the interval 


10,7} since. e{k) 1S symmetric with respect. to the pornt 


on = W.  ,lt wevlet, 


i = X= x (3-43a) 
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then x and y are to be determined from the 


(3-43b) 


equations, 


E(k.) = 0 (3-44a) 
aes) = 0 (3-44b) 
Or, TLOMm serting Ce =1 and using vequatvons) (3-37) fare 
(329). ; 
F, (x,y) cos“x + F, (x,y) cosx + BF, (x,y) = 0° (3-454) 
F, (x,y) cos*y + F, (x,y) cosy + By (x,y) = 0 (3-45b) 
where, 
Fy Cs,y) eit. Gc. 1) 
= -a{2(y-x) - 7 + (sin2y-sin2x)} (3-46a) 
Fo (x,y) = 1/2 - btn 2 OL) IR 
= is - 2(1-20) (siny - sinx) } (3-46b) 
F, (x,y) = (14+2a)S 
= (1+20){2(y-x) - 1} (3-20c) 
Eliminating the physically unacceptable solution x = y, 


i.e. I = 0, from the system of equations (3-45) we obtain, 


“te 


- ! : . : 
i a 


(tep—ey 


” a . 7 mar os Le 


om and 


vantolteapa oda ort 


(akb-E) 


(4hh~o4 


. : ; ; : 
‘ep a! ae. “7 WOLFSEUBS PILE? irises i= 7™ re “tis 
® 
- 


(vYixe).v + =z ODTXe) oF t «Root cat) oF 


+ Vv \~aXtia a : a ana (ys x) r@ 


jix- = (ys) 
| ae P 
rf c\e = ty ae ae 
) ae © te = : 
Blokel! <= (ty oe 
=e LS eee) 2 ter 


- : 
=o : 


o 


4 


.¥ * Aolgo.Gn slasiqsoosny vbiawtewia gid sg 


Logdo aw (cé-E} agoragps to medayn SAD aos? ,0. 8: Bee 
" ' 


ve ie 


F, (x,y) cos*x + F, (x,y) cosx + F, (x,y) = 0. (3-47a) 
F, (x,y) fcosy + cosx} + Fy (x,y) =) (3-47b) 


Define the following: 


a = oe (3-48a) 
bs i> (3-48b) 
Thus, 
y = atb | (3-49a) 
x = a-b (3-49b) 


Since x and y are restricted to the ranges 0<x,y<7, then 
from (3-48) we have O<asn and =1/2<b<7/2. In terms of 


the symmetrized variables a and b, system (3-47) becomes, 


n, (b)cos*a 4: n, (b) cosa = F/2o=" 0 a (3-502) 

a{m,(a,b) - 4b}m,(a,b) + m, (a,b) = 0 (3-50b) 
where, 

n,(b) = léasinb cos“b (3-5la) 


ny (b) = 2(1-20)sinb - 2a(1-4b)cosb 


- 8asinb ene°h (3-51b) 
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and, 


m, (a,b) = kcos(a-b){m - (1-20)4cosa sinb} 
eee) ia) (3p 2a) 
Bs 2 
m, (a,b) = cos” (a-b) (3=52b) 
m3 (a,b) = T =~ 2co0s2a SsinZzb (3-52c) 


The system of coupled transcendental equations (3-50) has 
been solved numertcally for various values of a in the 
range -w<q<-%. We have looked for solutions a and b in 
the ranges O<a<q and -7/2<b<7/2, that is for solutions 

x and y in the ranges Q<x,;y<i.) Knowing x andy we then 
Calculate the upper,bound forthe. exact,ground.state 


energy per site, Ey LY equation (3-26) 


Fgh [U-T+T?-U" +yt0 (V- THT? +2U7-Vv"-27TV 
e)} 2 ER (3-53) 
where, 
T = (y-x)/t (3-54a) 
U = y(siny = sinx) /1 (3-54b) 
Vee (ein2y = ssin2x) fon (3-54c) 
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For each value of a there are at least two solutions for 
Seconda im the interval 10,7.) Thewepucious solutions 

are eliminated by making use of the result that, since the 
one-particle energies, e(k), must be symmetric about dy = 
Terme 1nterval tf must also be synmetuicy about >} = 1. 

The results of oux calculations, which are identical with 


those of Niemeijer??, are given in Table l. 


Pa which is 


In Figure 10 the trial energy per site, E 
an upper bound for the exact ground state energy per site, 
Eyr LSeplotted as a function ef on Jy ars aAtso Gin 
cluded for comparison is the upper bound derived by 
Majumdar and Ghosh*/, namely By =o 7 lo) >) their idetinis 
tions of Jo and a are Slightly different from the ones 
used here. For a20 both approximations give the exact 
ground state energy per site, for -%<a<0 the same upper 
bound, but. for o<-% our upper bound gives a much better 
result. The first derivative of the trial energy per 
Site with respect. to. fis. continuoustar q = -—4> the 
second derivative, however, is discontinuous: We see 
that the nature of the upper bound for the exact ground 
state energy per site changes at the same value of a, 
namely a = -%, aS it does in the classical Heisenberg 
chain. dm this» approximation then, <q = °-% is the ipomt 
at which the antiferromagnetic second neighbour interac- 


tions become strong enough to alter the ferromagnetic 


ground state. 
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Table 1: ‘Numerical solutions of equations 


(31= 50) for varzeussvalues of o, 


f 


with corresponding Ey 
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Table l 


x (radians) ~y (radians) Ey 
0.215914 0.255793 -0.185005 
0.469102 0.654895 ~0.175500 
0.632849 Teg23o7 =Ol661 71 
0.695777 1.30736 ~0.161076 
0.719143 1.47114 0.159422 
0.736543 1.57080 -0.159818 
0.752498 Ticecon =0.161506 
0.767143 1.68966 -0.164093 
0.780414 1.72948 <0. 167351 
0.792365 1.76194 -0.171129 
0.803112 1.78912 ~0.175323 
0.812790 Lami OR 79657 
0.821530 Dees 246 0.184675 
0.829449 1.85013 EG. 189731 
0.836652 1.86580 -0,194991 
0.843227 1.87982 -0.200426 
0.886869 1.96790 -0.260779 
0.910043 peonoAd 0.326636 
0.934246 2.05799 0.464345 
0.954467 2.09564 -0.747501 
0.970220 ey, -1.46508 
0.984898 2.15196 -14.4573 
0.986394 2.15472 Sit aot 
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Figure 10: 


Jpper bounds for the exact ground 
state energy per site for an infinite 
chain with J,<0, aS tai (RUNCELOM Of a. 


For Niemeijer's result, asymptote 


for large negative a is plotted. 
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We now want to find the asymptote to the By versus 


oC curve LOT J, <0 (Figure 10) \fer) Targe necative a: First 


we define, 
Liens 
eye a | (3555) 


Moin cnerract. thaat < = s(6:)); -ye= Vip ape ane dividing 


(3-47) by oAawe. optains 


Fle (eee) poe 0 (3-56a) 

g(x(8),y (8) ,B) =) (3-56b) 
where, 

E(x,y,8) = G, (x,y) cos*x + G, (x,y) cosx 

+ G, (x,y) (3-5 7a) 

G(xe¥ 76) = G, (x,y) {cosx + cosy} + Gy (x,y) (3-57b) 
and, 

G, (x,y) = BF, (x,y) iia (mis Ne) . (3-58) 
Define, 


BE, 


BE = E(x(8),y(8) -8) £ 


(3-59) 
where Ey is given ineequation y(Ss53)inctwhen @i-fsowelL.e. 
gesn07, xv’and ‘y (also a and b) approach certain values, 


Seven ond, ye ee ran De y. “OMultiplying through» équations 
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(2-00) by "sand Vetting’ '— +7 0" wemhave, 


(1 - 2cos2a,sin2b, - 4b,)cos* (a,-b,) 

+ 4cos (a -b ) cosa tsinb re Gent Seu OES eA Co odent he 
16sinb_cos“b, cos?a_ eho | 72:(n-4b.) cosh) 

+ 8sinb,cos“b, }cosa, SO (3-60b) 


From these we obtain, 
cosas =O (3-61a) 


4sin’b, ~ (7=4D.) cosh =F 40) (3-61b) 


whose solutions in the allowed ranges of a and b are 
oer w/2 radians, bs = 0.534236ered1ens., Phen from 


(347), x, = 0-38 po60 radians and yi ge) 2.15506) radians . 


We expand x(8), y(8) and 3 rad ese Wake aa (che Late a_Taylor 


Series about 6 =) 0" as toblows: 
a Baier e ahead A,B + 0 (87) eae ae q (65.624) 
Vite Ven oe B,8 te 0(87) Her iran (3 =6.25) 


B(x (B) ,y(B) ,B)= Cy + C,8 + Oe le Ps eo) 


0 


Not knowing = anc aigexpiicitly as functions: of By, we 


cannot find the coefficients Ay and B. by direct. dififer- 


a 


entiation; instead, we use equations (3-56). Substituting 
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a” aon chibies 
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io} « pinky rink 
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. aces jiel 1rt® tte 


the first two of equations (3-62) sinto (3-56), we expand 
the left hand sides of the two equations in a Taylor 
Serice ah 6 about 6 = 0 .- Then, equating coefficients of 


equal powers Of 8 ‘on both sides, we, find ‘to First ‘order, 


ae Jf of st y 

E Ay aa See + 7 we: =) (3-63a) 
B=0 

9g 5 cs) ols) = 2 

[Be a, ao ay By +: 5 ae 0 (3-63b) 


Solving this system of two linear equations in the two 


unknowns A. and B. we have, 


a ie 
FF] seo- (98) ano- 7 (BE) peo (8H 
SRE ee ed (eo s Aare Neen cola 2 

p= ee el eGo) 

Ca as 
peg Poe Te RS ine es Baie 
FF) sao (3) gaom 7 (Fe peor lie] 

x 2B) p97 19%) g=07 9X) g=97 198) g=07 (3-64b) 


ti we also substitute the first two ofvequations (3-62) 
into the third and expand the left hand side in a Taylor 
series in 8 about @ = 07 and compare with the right hand 


Side, Elen. te 4firset order, 


coe re al © ter ger RD, (3-65) 


and, 
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Wer Can. rewrite, (3-62c),aianrthe form, 


Saleen 9 ues al Ge 


i: 0 it ON e oe tate (3-677) 


From (3-67) we see that the equation of the asymptote for 


large negative a is, 


t = — 4 
Ey = acy + Cy (3-68) 


sO Co is the slope and Cy is the Ee intercept of the 


asymptote. Performing the numerical calculations we get, 


ij 


x (8) 


DA 55 Ode OMS Ene ten Oe ie eet o 26on) 


y (8) 


E(x(8),y(B),B) = 0.144405 - 0.01640428 


PROKIC MV MEE ih aWieaGOc) 
and the equation of the asymptote is, 


By = 0.144405a - 0.0164042 Go 708) 


The asymptote is shown in Figure 10. 
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B. Antiferromagnetic NN Interactions 


For J5>0, nearest neighbour interactions are anti-— 


ferromagnetic so we take Jy = +1 and therefore a = +55 - 
We want the solution of (3-39) for Jy = +1, -“<a<+o, 
Niemeijer*? has Said that I =aa/2,1jeus the solution for 


all valueseotia. This, however, ,16esnot correct, . 
Tf I = {n/2,7r} is the solution then from (3-37) we 


havess =“Qe="0). R-="-2." Equation’ (3-39) becomes, 
e(K)wse (p+ 2 (1-20)) cos¢ (3-71) 
T : K 


Bom jastn74~ (2i. 28540) "the coer ivcient (or "COsp, " nS 7p) 
is positive, for a>s%+n/4 it is negative and for q = %+7/4 


RewtcuZere:. 


Thus tae fT /21s ase: SOluewon or (3539 itor o<s+1/4 
when J, = +1 but, for C2sti/47) lean len Cannotube 4 


SQUlUtLLON Since, if 1t is, then) eth) 2.0 when 6, C1 and 

e (k)<0 when dE T, which is contrary to what we assumed 
about the properties of the solution interval I. For 
a>st+n/4 we look for a different solution interval. 
Niemeijer's results have that the upper bound for the 
exact ground state energy per site approaches +o as a > 
+o, Since the upper bound approached -~ as a > +© for 
the ferromagnetic nn case, Stephenson~” suggested that 
Ey should exhibit a similar behaviour for the antiferro- 


magnetic nn case. This speculation led to the discovery 


of the error in Niemeijer's results. 
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(3) J. = FIM soe 007 0} dg dO Seti as (mee yRevis the, sol= 


4 
wtzon. Thus, SG = Q = 0, R= =-2. The elementary excita-= 
tion spectrum (single-particle energy) is given by 
equation (3-71). The upper bound for the exact ground 
state senergy ‘per ‘*sitesis, from (3526), 


1 a 
eal t(,- ' | 2 
‘ eae’ = 20) |. 2 Eo | (S72) 


So far we agree with Niemeijer's calculations. 


UA) he = +L, Daoete 20, J5>0) : For 0<a<44+71/4 a possible 


1 
POlLULLOn 1S Stl) given by P= in/2)1 in whieh cace 
equations (3-71) and (3-72) give the elementary excita- 
tion spectrum and trial energy per site, respectively. 

For o2%4+1/4 the interval I is now {x,y}, as shown in 
Figure.11. Again we restrict ourselves to the interval 
{0,7} since e({k) is still symmetric about Oper Follow- 


ing the same procedure as in Case (2) we obtain, after 


eliminating the solution x = y, 


F, (x,y) cos*x + Fy (x,y) cosx + F, (x,y) = 0 (3-732a) 
F, (x,y) {cosy Hrcoss pas Fo (x,y) = 0 (3-73b) 


i Fos F, are defined as before. System (3-73) is 


identical with system (3-47). If we solve (3-73) numer- 


where F 


Leal ly for x -aud yy we can then find the trial energy per 


site from (3-26): 
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where TIT, U, V«are defined by (3-54). 
System (3-73) has no solutions in the allowed ranges 
Oras ang Vj, 02s, Von, LOrv- “sos but does, have valid sol 
Utlons for all a@Z1. We have already Shown that I= 
ie on) LS: a= DPOSsslLpLessolLution Of equation (3-59) Lor 
O<a< %+7/4. The actual solutions of (3-39) for Case (4) 


thus are: 
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{x,y}, where x and y are determined by the 
system of equations (3-47), for a>s+1/4. 
cc). Huther-i=%i/2ontoun Bee 2 iy; ite weer has the 
Tower traialweneroyspem site, for 120<447/4. 
Ey Eom oa 2) one Et trom (3=/74) intersect when aq 2) 1.05. 
For 1$a$1.05 equation (3-72) gives the lower trial energy 
Ber cite while for 1. 05s0c+1/ 4 deere (3-74). “Thue, i = 
i,f2o7, ag the solution intervaleter J20.1205 and tor 
Wee o<eta 7 ote ser Silk. 
Numerical solutions of system (3-73) for various 


values of a in the range 1.05Sa<+@ are given in Table 2. 
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Figure 1: 


Elementary excitation spectrum, e(k), 
asfrunction oF by for different ranges 


of a with J,>0. 
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Table 2; Numerical solutions of equations 
(3 - 73) for various values of o 


with corresponding Ey 


76 


Cooker 
ory Lay Ay Wak (EY = &) 


nin eye 1 IoD rH S| iw 


1000 


ma 


Sera. 


a2 5 OW 


de 1SPO33 
veto Ox. 


Lease rel) 
Mel PSiyesie 
1.12401 
dee Oe 
Leos 2 4 
Veo 523 
LL08826 
By tres g7An 
1.06043 
1.04686 
1203074 
1.02144 
1.00944 
PeOn3ae9 
0.98823 


Oe95 Gy Z 


Table 2 


ans) 


os 


uy 


7 


y (radians) 


2.58948 


ZNO d L 
2etd eo2 


Bag 7? 
2.43884 
o041233 
2.39080 
Pesjesy 
2.35764 
2.34451 
2.33306 
229217 
2.26682 
2.23682 
2.21957 
2.19734 
Behe ool 
Peeode 


an U5 534 


= 


=-0 


=) 


t 
Be 


e ALIS S 


HaOGH 1 
Se Ala 


LOO he 
PALO SUG 
eo) 
.250458 
weOZoG2 
(2a LG 
o20 7007 
AAS ho WheKs: 
Oo ou 
~-434077 
so OCT 
~ 715597 
SOA ES al By 


A3¢04 


“14.4245 


~LAe 380 


ei 


_ 
- fA 
ne 
: 
ify 
7 | 
ul 
7 7 


a (anak Bas) vo 


a 


RE eCOS  H- BROSe.S 


(Taeos  0= - WTsres 


em : a ; : 

Teaae.f 

; as’ ae 

os 
a " a oe 

O£BS9)..0 te ge 9a0Edt | ae 

0! i CERTE.S £0stiet ~ 


protre. b= CAnth.o 


BIS. 0 cvs \s.s 


(in HS80GF. Ss : (VSLived 
TseoLl 


0 reese ,< ASHOO. 2» 
ft. Slant 
cises.¢ conan 
ch a gaa een. 
cs.3 BTCEO LL 
firs s PRICED. SE 
evr. 4 BBOOD. TF 
f2aRE.S eeeoa.f 
bSh hI] feaics fESIRE OD 
~ bal bEczL.S Vaveee .o 


In Figure 12 we have plotted the trial energy per 
site, shor aS vaaetunction of cy; Jy = +1, for the results 
of Majumdar and Ghosh~/, Niemeijer?? and ourselves 
feaquatilon (6-72) 4and Table 2)... "Ferva = 1.05 we have 
rederived Niemeijer's upper bound; however, for a 2 1.05 
our upper bound is considerably vera The upper bound 
for the exact ground state energy per site that we ob- 
tained as a correction to the work of Niemeijer is con- 
tanuous, at the point of antersection, ore l.05; but its 
first derivative with respect to a exhibits a disconti- 
nuity there. Niemeijer's and our upper bounds are both 


much better than that given by Majumdar and Ghosh. For 


QO: (037 2M. guves Ey (a=0) =O) 419630 ey and tines 


0 
agrees with the result of Bulaevskii>’. The well known 
exact result due to Hulthén-- is Eg (a=0) So ein hy 
e -0.44315. 


Finally we want to find the asymptote to our By 


versus,.a curve for J, >0 (Figure 12) .for large positive aq. 
We proceed as we did in finding the asymptote for the 


case J, <0. Gompar ence o( 3-55 )s wanda 4 Je, 


(EL) ts BG) (3-75) 


WREN .0.0% chop Lei a.8, 420 eand wy g@pproach and yy, 
Pespectiveiy, where x = 0,966S5200raddans and y= 
2.15506 radians, as before. We make Taylor expansions of 


(6), oy Ch). and E(x(8),y(B),8)- Proceeding as for J, <0 
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22 | en ee EE | a 8 (37763 
i’ B=0 iy p= Oi 
J,>0 J, <0 
where x) Srey Xo = V5 X2 = By and, 
si ecb act ha ie 
(on att ae Eee re On ere tm ore) 
4 a: 
We find, 
on) 50 Spine) lige-o ake) 
il ié 
ae sda: 0 (3-78b) 
and, 
SSW Riea arc Fee ae eats) 


The Taylor expansions for x(f8) and y(8) are identical 
with (3-69a) and (3-69b), respectively. However, for 


E(x(8),y(8),8) we have, 


E(x(8),y(8),8) = -0.144405 + 0.01640428 


Bo a) fe 0 a (3-80) 
The equation of the asymptote is found to be, 


Ey = -0.144405a0 + 0.0164042 (3 ils) 
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This is just the negative of the asymptote we had for 


J,<0, equation (3-70). The asymptote (3-81) is shown in 
Figure 22% 
Note Added In Proof: For the ferromagnetic nn case I = 0 


is the solution for a>-k. By inspection we see that I = 
0 is a valid solution interval foe large ferromagnetic 

second neighbour interactions when first neighbour inter- 
actions are small, and of either sign. For large ferro- 


Magnetic J nnn spins will be aligned. In thes present 


Zt 


approximation, nn spins are aligned also, even when J,70. 


For large ferromagnetic Jo the excitation spectrum is, 
= 4 1 = 
E(k) = 2a (cos$, 1) (cos, t1+5=) (3-82) 


and the trial energy per site is, 


(oa ak ' = 
Ei +k (lta) > Eo (3-83) 


Bavatzons (3572) vanes (3-355) 2£0n Ey intersect when a = 


-14,,14. 5B, un (3-83) has the correct asymptotic slope as 


t 
i 
a> -o and is included in Figure 12. 
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Chapter VCore lusion 


In this thesis we have introduced and studied cer- 
tain properties of two one-dimensional lattice models. 
For the spin % one-dimensional Ising model having both 
nn and nnn interactions in zero magnetic field, the canon- 
ical partition function has been found using the transfer 
matrix method and an equation for its locus of zeros was 
Then OoLal Meo tas, a function, ob MM, the 2ario of the next 
nearest to the neerest neighbour interaction energy. The 
equation was solved numerically for different i values 
and the results were plotted in Figures 1-6. We then 
discussed general properties of the loci of zeros and 
investigated exactly soluble cases. 

Also considered was the one-dimensional spin % zero 
field Heisenberg model. Reviewing the work of Niemeijer, 
by means.of a temperature dependent Hartree-rock approx- 
imation we were able to derive an upper bound for the ex- 
act’ ground state energy per site at zero temperature for 
various different values of the nearest and next nearest 
neighbour interaction energies. After Seed that part 
of the upper bound given by Niemeijer for the antiferro- 
magnetic nearest neighbour case (J,>0) was wrong, we were 
able to obtain the correct upper bound, which was consid- 
erably better than that of Niemeijer. With a = Jo/Iy1 we 
also found the asymptote to the upper bound as a function 


of a for large negative a with J,<0 and for large positive 
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o® with J,>0. 

As indicated by Niemeijer this variational method 
that was applied to the zero field Heisenberg model can 
be quite readily extended to the problem where there is 
a non-zero magnetic field present in the z (Cartesian 
coordinate z) axis. 

In this thesis we have introduced next nearest 
neighbour interactions into two one-dimensional lattice 
medelshintansattempt” to bring *thém perkaps tas livele bit 
closer to describing an actual physical system. The 
solutions we have obtained will, it is hoped, shed 
more light on ideas about lattice systems and on as yet 
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Appendix 


In this appendix we discuss boundary conditions and 
the locus of zeros of the. partition function (2-4) in 
Chapter 2 in order remove any ambiguities that may exist. 
Equation (2-4) gives the canonical partition function of 
an N spin % one-dimensional Ising linear chain having 


both nearest and next nearest neighbour interactions pre- 


sent in zero magnetic field. We rewrite this as, 
Lc a 1 N-1 
Za (Ky 1K570) = a ee exp (K, } Cuca 
Oy=-l OWE i=1 
iN Sa 
+ Ko ) 0.9, ,5) (A-1) 
i=1 


where the superscript "fc" indicates that linear chain 
boundary conditions are used. Dobson's result, (2-5), 


becomes, 


gC 1,42 
Zy (Ky 1Ky/0) = 22,7, (Ky /0,K 


: ) (A~2) 


if 


Mo evaluate.the N-1 spin partition function, (K,,0,K)), 


ay-1 
we used the transfer matrix approach; however, it is well 
known that to use a transfer matrix in the manner that 
we did'‘in equation (2-11), one. must assume cyclic boun- 
dary conditions, i.e. the linear chain has now become’a 


ring. We then set this N-l spin function equal to zero 


to obtain its locus of zeros so we actually found the 
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(K,0,K,) = 0 (A-3) 


woeerve thes superscript "ring”™ means, that cyclie boundary 
conditions were employed. 


We now derive a relation between the partition func- 


Gion Of a ilainear chain and’ that vor a ring. | We have, 
Pees 0K sy VL eee en (A-4a) 
N-L* 28 ui cas ii 
all 
states 
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We notice that, 
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all 
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(A-6) 
Since exp (- Ko ey On)= coshk,, - On oysinhk, this becomes, 
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‘ ring ( 
~ SinhKy Oy 1% y-1 (Kp19+Ky) J 
(Aap 
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where cee on Peers ace K ie is the two-spin correlation 


function between spins at sites N-l and N (i.e. site 1) 
on a ring Of N-1 spins. Because of translational symm- 
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ring Cing vapials fy) Z 
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w(K, ,0,K ) = ©1955 ce Che Oe (A-9) 


(A=-7) can be written as, 
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(K,,0,K,) tcoshk, 


- sinhKjw, (K,0,K, }} (A-10) 
Now, 
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xX = Magnetization = 1-4 (Aaa) 
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Equation (A-10) thus becomes, 
- sinhk, (x“+y“y) } 


Ce. ae _ .N-1 
Zo (Ko 191 Ky) = ma {coshk,, 


WN-1 tcoshk = sinhk, (x“+y/y) } (A-15) 


ves 2 


Setting (A-15) to zero we obtain as our locus of zeros, 
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at ct (-1) 17N-11/N-1 (A-16) 
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where, 


coshKk, - sinhK, (x*+y*/y) 
D = i RR LE (A-17) 
coshk, - sinhk, (x t+y y) 


Recall’, 
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We can rewrite (A-16) as, 


A : 

ae ayer 0) yee *exp| dargp| (A-19) 
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We are “interested in the case Neo. Fortthis,. (=1) 
becomes continuous along the unit circle, fay ee Gea 
and exp|fargp] +L. VThus); as Ne 

ae, 1/N-1 

La en (A-20) 
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(A-20) is the” locus of zeros,-as "N’>+*o, obtained’ from, 
Ge om a 
Zu (Ko 19K) = 0 (A-21) 


The equation that we actually solved numerically was, 


EVOMent 2). 


Ail 


a Taira Ja! - 


) iA Sd | gh 
£ Linen ing! UNE |g APs [ 


—— 


we - a 
(4 +r" , ~ |W ge@an oat ns hnatseoyotnt S36 5 J 
~ a. 


2 7 7” 4 ‘ ——. 
tev sit paols. apegaisnes Bae 2O9¢ 
, 

7 


aB our? 9 i aad laexag| qn ; ts 


» 


1 


= 
A 
~“\ I 
i~| : (t=) a ie 
} i 
al r 
tisjdo i @6),e01ss Io aesok eae es (OSA) 


iyenut bovier' ¢iisugse oe dens ees 


wiLech-ls the vexact “locus “of zéros"f6r"all N “obtained from, 


ring 


24 kan 


Rie ny Ke 


Xo eaye (Aes) 


al 


DaAus we see Chat ounm,procedure of getting the locus of 


ve : ue Ke 
N (K, 7K, 70) = 2a 41 


by solving for the zeros of the ring (z (K,,0,K,))} was 


zeros of the linear chain (2 Cha) 
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